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Abstract. We have experimentally observed re-equilibration of a magnetically 
trapped cloud of metastable neon atoms after it was put in a non-equilibrium state. 
Using numerical simulations we show that anharmonic mixing, equilibration due to 
the collisionless dynamics of atoms in a magnetic trap, is the dominant process in this 
equilibration. We determine the dependence of its time on trap parameters and atom 
temperature. Furthermore we observe in the simulations a resonant energy exchange 
between the radial and axial trap dimensions at a ratio of trap frequencies LUr/cUz = 3/2. 
This resonance is explained by a simple oscillator model. 



PACS numbers: 05.45.Xt, 39.25.-|-k, 05.30.Jp 
1. Introduction 

Magnetic traps have become a standard and inexpensive tool in cold atom physics in 
recent years ff . Together with the application of cooling techniques such as evaporative 
cooling and sympathetic cooling, they have enabled the formation of Bose-Einstein 
condensates and, more recently, degenerate Fermi gases 0. The most commonly 
used species for cold atom experiments (Rb, Cs, Na) are easily cooled to the micro- 
Kelvin regime, where the atomic dynamics in a magnetic trap is usually well described 
by regular harmonic motion. For other species such as the noble gases jH |3] and the 
group II elements [HI EI however, the lowest attainable temperatures can be in the milli- 
Kelvin regime. In this regime the atomic dynamics is more complicated, with coupling 
between the different dimensions of the trap becoming important. 

Anharmonic mixing is a process that couples the motion of atoms in different 
dimensions of a magnetic trap. This coupling enables redistribution of energy over 
the coupled dimensions. Therefore, anharmonic mixing has to be considered when 
applying techniques that rely on the motion of atoms in a trap, such as evaporative 
cooling |H], Doppler cooling in a magnetic trap jH] and rethermalization experiments 
aimed at measuring a scattering length [TU]. In order to investigate on what timescale 
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anharmonic mixing plays a role, we perform a Monte Carlo simulation of a cloud of 
atoms in a loffe-Quadrupole magnetic trap. We determine the dependence of this 
timescale on trap parameters and on atom temperature. Surkov et. al. ^T] have 
investigated anharmonic mixing in the limit of low temperature where it can be treated 
as a perturbation of regular harmonic motion, but to the best of our knowledge this is the 
first paper discussing anharmonic mixing at higher temperatures. We will also present 
here an experiment showing that anharmonic mixing does play an important role in the 
atomic dynamics, and can even be the dominant mechanism of energy redistribution 
under certain conditions. To conclude we will look at possible applications of anharmonic 
mixing. 

2. Monte Carlo calculations 

A loffe-Quadrupole magnetic trap (MT) is one of the most common magnetic traps. Its 
potential is given by 



where x and y are the radial trap dimensions, z is the axial trap dimension, /x is the 
magnetic moment of the atom, a the gradient of the magnetic field, (3 the curvature of 
the magnetic field and Bq the magnetic bias field. In this expression for the trap potential 
terms of order higher than four have been neglected. The MT has trap frequencies Ur 
and uJz in the limit SksT <^ /i-Bo^ where the trap shape is harmonic: 



Here m is the atomic mass, r = (x'^+y'^Y^'^ the radial coordinate, cUr = [■^{a'^ / Bq— P)]^/"^ 
the radial trap frequency, and Uz = (yU/S/m)^/^ the axial trap frequency. The higher 
order terms in (P) are in this limit negligible compared to the harmonic terms, the term 
2a(3{x'^ — y'^)z that couples the motion in the axial and radial directions is therefore 
absent from (j21). 

The starting point of a simulation is an atom cloud in equilibrium in a harmonic 
trap. The initial positions and velocities of the atoms are determined by a Monte Carlo 
method. For the properties of the atom we use the values for metastable neon in the 

state, as this is the atom we use in our experiment. Two kinds of clouds have been 
used, namely thermal clouds and clouds in which all the atoms have the same energy. 
The former is used to compare the results of the simulations with experiments while 
the latter is used to gain a better understanding of the dynamics of atoms at a certain 
energy. During the simulation the position and velocity of each atom is determined as a 
function of time by integrating the equations of motion. At the start of the simulation 
the magnetic bias field Bq is adiabatically ramped down, thereby increasing the energy 
of the atoms in the radial direction and changing the effective shape of the potential 
from harmonic to almost linear. The ramping is adiabatic if the condition < 1 



U{x, y, z) = fi[Bl + {a^ - B^l3){x' + y^) + 2B^l3z^ 

+ \p\x' + y^f + + 2aP{x^ - y')zY'^ - ^,B 



(1) 




(2) 
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Figure 1. Monte Carlo calculation of the radial temperature as a Junction of time 
for a thermal cloud with an initial temperature of 1 mK and lO"* atoms (gray). The 
exponential fit is also shown (black); the temperature during the ramping of the bias 
field (t — Q — ms) is left out. 

is fulfilled jTTj. The simulation yields the kinetic and potential energies of the cloud in 
every direction as a function of time. From the kinetic energies in different directions the 
transfer of energy from one direction to the other can be determined, and a comparison 
with the experiment can be made. 

We choose the trap parameters at the start of the simulation as a = 1 ■ 10^ G/m, 
(3 = 47.5 ■ 10^ G/m^ and Bq = 99.6 G. The bias field is then ramped down to 1.5 G. 
These parameters are chosen because they are easily accessible in our experiment. A 
typical calculated result of the temperature evolution, or the average kinetic energy, in 
the radial direction of a thermal cloud of 1 mK after ramping the bias field is shown in 
Figure H The first few tens of milliseconds show a linear decrease of temperature, after 
that the decrease is approximately exponential with a mixing time of 112 ms for this 
particular cloud with a temperature of 1 mK and a relative change in temperature of 
8%. The linear decrease of temperature in the first tens of ms is caused by high-energy 
atoms that mix on a timescale of ^uj^^ = 11 ms. This is the average time an atom needs 
to reach a singular point in the potential ^Hj, where almost instantaneous mixing can 
occur. 

The dependence of the behavior of a thermal cloud on initial trap parameters and 
cloud temperature is relatively weak. When a is varied over the range 0.9 — 1.2- 10^ G/m 
and Bo^end over the range 1.5 — 20 G, the mixing time stays within the range 100 — 130 
ms. Only when (3 is decreased, a significant change in mixing time occurs: Tmix ~ 400 
ms when (3 = 37.5 ■ 10^ G/m^. The reason for this weak dependence on trap parameters 
is that a significant part of the energy of the cloud is carried by atoms that are in a 
regime where mixing does not occur or where it occurs on a timescale of \^^^ ■ 

To show that indeed the timescale of energy transfer between dimensions of an 
atom depends strongly on the energy of that atom, we now take clouds of atoms that 
all have the same energy Eatom in order to determine the mixing timescale as a function 
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Figure 2. Calculated mixing time for atoms with a fixed energy Eatom drawn as a 
function of Eatom/Sks- 

of energy. This gives us the average behavior of an atom with a certain energy. The 
result is shown in FigureEl Below Eatom/^kB = 375 fiK no mixing occurs on the longest 
timescale we considered, i.e. the lifetime of the metastable state Ne*(^P2) of 14.7 s |12j . 
For temperatures above 1.5 mK the mixing time is determined by ^uj^^ and becomes 
independent of energy. 

Next, we investigate the dependence of Tmix on the trap parameters a and /3 
while keeping the energy constant at Eatom/^kB = 700 /xK. Figure El shows that Tmix 
increases with a and decreases with p. This behavior can be understood qualitatively 
by comparing the strength of the harmonic (~ — f3) and coupling (~ af3) terms in 
©• 

Furthermore we determined the mixing time as a function of the value of the bias 
field at the end of the ramp. When the final value of the bias field is higher, less energy 
is added to the atoms and the final temperature of the cloud is lower. Therefore we 
expect that the mixing time becomes larger for higher final values of the bias field. As 
can be seen from the result in Figure ID however, there is a broad resonance (meaning 
small mixing time) in the mixing time around 14 G. This occurs exactly at a ratio of 
radial and axial trap frequencies of ujr/^z = 3/2. The resonance is so broad because 
the radial trap frequency changes only very slowly with bias field. The reason for this 
slow change is that at higher bias field the final energy of the cloud is lower, resulting 
in a higher oscillation frequency in the MT because the trap is not harmonic, but at the 
same time the trap is less tightly confining, resulting in a lower oscillation frequency. 
These two effects almost completely cancel. 

3. Resonance model 

In order to explain this resonance the dynamics of an atom in the MT is now described 
as two coupled one-dimensional harmonic oscillators, whose energies represent the radial 
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Figure 3. a) Calculated mixing time as a Junction of magnetic gradient in radial 
direction a. h) Mixing time as a function of magnetic curvature in axial direction (3. 
Atomic energy is Eatom/Sks = 700 fiK. 
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Figure 4. Calculated mixing time as a function of bias field after ramping. Atomic 
energy is Eatom/^kB = 700 jj^K. 



and axial energies of the atom. The couphng term is chosen to he W = x'^z, as shown 
in Figure Now we examine the forces and that the two oscillators exert on one 
another, assuming they oscillate as a; ~ cos^ujxt) and z ~ cos{uJzt): 



Fz = --^W = x^ r-u cos^{ujJ) ~ cos(2u;^t) + 1, 
d 

Fx = = ~ COs(co'a;t) COs(ti;2t) ~ COs{{Ux + OJz)t) + COs((ci;2: — LJz)t). 



(3) 



A resonance occurs if an oscillator is driven by a periodic force with its 



Anharmonic mixing in a magnetic trap 



6 




Figure 5. Schematic picture of the coupled oscillator model. Two one- dimensional 
oscillators with frequencies Ux and lOz can exchange energy through a coupling. 

eigenfrequency, Fj ~ cos{ojit). From Q it can be seen that this happens if uJz = ^u^- 
In general, a couphng term x°'z'' yields a resonance at au^ = hu^. When potential ((H) 
is linearized a term x'^z^ appears, explaining the resonance in Figure |3] A more formal 
treatment explaining the resonance is possible by transforming the Hamiltonian of an 
atom in the potential to Birkhoff-Gustavson normal form but that is beyond 

the scope of this paper. 

4. Experiment 

To check our simulation results, we now compare with experiment. The only difference 
between simulation and experiment is that in the experiment the atoms can also collide, 
which we did not include in the simulation. Therefore it is necessary to consider both 
anharmonic mixing and coUisional equilibration of energy. Our setup has been described 
previously in |S1[T3]. Briefly, a discharge source is used to create metastable neon 
atoms, and after passing through several laser-cooling sections the atomic beam flux is 
6 ■ 10^'^ s^^ [5]. From this atomic beam about 10^ atoms are loaded into a magneto- 
optical trap (MOT) at a temperature of 1 mK. These atoms can be transferred to 
a loffe-Quadrupole magnetic trap with almost unity efficiency. To match the size 
and temperature of our MOT the trap parameters are chosen as a = 1 ■ 10^ G/m, 
13 = 47.5 ■ 10^ G/m^ and Bo = 99.6 G. The MT has trap frequencies of Ur = (27r)376 
Hz and = (27r)44.7 Hz. 

For this experiment the MOT was operated at a detuning of — 2r and at an intensity 
of 0.5 mW/cm^ for the three MOT beams together. After turning off the MOT, a 50 
fis long spin polarization pulse of 2 mW/cm^ was given to put all the atoms in 
the \mj = +2 > state, and then the MT was turned on. The atom cloud now has a 
temperature of approximately 1 mK and contains approximately 1 ■ 10^ atoms. After 
turning on the MT the atoms are held there long enough to allow the atom cloud to 
reach its equilibrium state. Then the bias field is ramped from 99.6 G to 1.5 G, adding 
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Figure 6. Radial rms radius of the atom cloud as a function of time. The initial 
temperature in the MT is 1.0 ± 0.1 mK. 



energy to the atoms in the radial direction and compressing them spatially. The cloud 
radius in the radial direction as a function of time after the ramping is determined by 
absorption imaging. The radius is directly proportional to the potential energy and 
therefore the temperature because the trap shape is linear. A typical result of a series 
of measurements is shown in Figure IHl 

The exponential fit gives a characteristic equilibration time of 130 ± 25 ms and 
a relative change in width of 5 ± 3%. This result agrees with the timescale of the 
simulation as found from the data in Figure H If the equilibration would be due solely 
to collisions, this would give a coUisional cross-section of a = 87r-(200ao)^, with qq the 
Bohr radius. However, the cross-section at 1 mK can not be larger than the unitary 
limit a = 87r(AdB/27r)^ = 87r-(80ao)^. Here XdB is the thermal de Broglie wavelength 
of the atom. Also, we did not observe a dependence of the equilibration time on the 
atomic density, indicating that anharmonic mixing is the dominant mechanism here. 

5. Applications 

An interesting application for anharmonic mixing is to use it to enhance the efficiency of 
Doppler cooling in a magnetic trap. This cooling technique is one-dimensional because 
only along the axial direction of an MT all atoms are polarized the same way. Cooling 
in the radial directions can be achieved by reabsorption of scattered photons JH] if the 
cloud is optically dense. If that is not the case anharmonic mixing can cool the radial 
directions, as shown in an early experiment by Helmerson et. al. [H]. By making the 
mixing time as short as possible the cooling can be fast and the atom losses as a result 
of the cooling can be limited. 

Another technique that could benefit from anharmonic mixing is evaporative 
cooling. One of the problems that can occur when an atom cloud is cooled evaporatively 
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is gravitational sag, in which gravity shifts the equipotential surfaces of atoms in a 
trap in such a way that they do not coincide with the surfaces of constant magnetic 
field. This causes the cooling process to become one or two-dimensional, reducing its 
efficiency jTTj. If the dimensions are coupled by anharmonic mixing the evaporation 
remains three-dimensional even when gravity plays a role. 

A situation where anharmonic mixing is undesirable and needs to be suppressed 
is in rethermalization experiments to measure the scattering length ^OJ QHl- The 
scattering length is one of the properties of an atom that determine the feasibility of 
evaporative cooling and reaching the transition point of Bose-Einstein condensation 
for a given number of atoms at a certain density and temperature [H]. It can be 
measured by observing coUisional equilibration after ramping the bias field as we did 
in our simulations and experiment. We can conclude that this method will only yield 
reliable results if the timescale on which anharmonic mixing occurs is long compared to 
the collision time. 

6. Conclusion 

We determined the timescale on which anharmonic mixing occurs and the dependence 
of that timescale on trap parameters and atom temperature. We observed a resonance 
in the mixing time as a function of the end bias field, and explained this with a simple 
oscillator model. We verified experimentally that anharmonic mixing does indeed occur, 
and that its timescale can be short compared to the timescale needed for coUisional 
equilibration. The application of anharmonic mixing to improve Doppler cooling and 
lower-dimensional evaporative cooling is possible. 
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